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We study a matrix product representation of the Bethe ansatz state for the Lieb-Linger model 
describing the one-dimensional Bose gas with delta-function interaction. We first construct 
eigenstates of the discretized model in the form of matrix product states using the algebraic 
Bethe ansatz. Continuous matrix product states are then exactly obtained in the continuum 
limit with a finite number of particles. The factorizing F-matrices in the lattice model are 
indispensable for the continuous matrix product states and lead to a marked reduction from 
the original bosonic system with infinite degrees of freedom to the five-vertex model. 
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Finding an optimal representation of many-body quan- 
tum states is the key issue for variational calculations. 
In the density matrix renormalization group (DMRG) 
method, 1 ^ which is a powerful numerical method for 
one-dimensional strongly correlated systems, the matrix 
product state (MPS) is used as a variational state. 2 ) 
In some cases, such as the Affleck-Kennedy-Lieb-Tasaki 
(AKLT) model, 3 ) it is possible to obtain the exact ground 
state as the MPS. While the MPS has a long history, 4 ' 5 ) 
it is currently attracting interest in the interdisciplinary 
field of condensed matter physics and quantum informa- 
tion science. 6 ) 

At hrst glance, it seems useless to find an optimal rep- 
resentation for eigenstates of quantum integrable models 
since they are exactly obtained from the Bethe ansatz 
wave functions. However, direct calculation using the ex- 
act eigenstates is a formidable task. A recent develop- 
ment in this field is the application of the factorizing 
.F-matrix (Drinfel'd twist), 7 ) which is a similarity trans- 
formation to a new basis (F-basis) in which actions of 
operators are simple. This new scheme enables a direct 
computation of form factors. 8 ) The Bethe ansatz has a 
long and rich history beginning with Bethe's solution 
of the spin-i Heisenberg chain, 9 ) which is now called 
the coordinate Bethe ansatz, of which there are several 
variants such as the algebraic Bethe ansatz (or quan- 
tum inverse scattering method) 10 ) and matrix product 
ansatz. The matrix product ansatz proposed by Alcaraz 
and Lazo 11 )gives exact eigenstates expressed as the MPS. 
Moreover, they have found algebraic relations between 
matrices constituting MPS, which reproduce the correct 
results obtained by the coordinate Bethe ansatz. That is, 
permutations of the spins in the coordinate Bethe ansatz 
are expressed by commutation relations of the matrices. 

In this sense, an MPS-type representation of the Bethe 
states is interesting in many fields of physics ranging from 
mathematical physics to quantum information theory. 
Motivated by this, we have studied the connection be- 
tween the algebraic Bethe ansatz and the matrix product 
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ansatz and have shown their equivalence in the Heisen- 
berg chain. 12 ) Surprisingly, the connection turns out to 
be related to the F-matrices. We have found that an 
MPS-type representation of the Bethe state itself is easily 
obtained by changing the order of the product of the L- 
operators. The concept is similar to the quantum transfer 
matrix 13 ) and will be explained in detail in the present 
letter. To find algebraic relations among matrices ap- 
pearing in MPS, the obtained MPS-type representation 
is too complicated. Thus, the F-matrices are necessarily 
required. In the new basis introduced by the F-matrices, 
the matrices have a simple structure. 12 ) This simplifica- 
tion corresponds to a mapping from the six- vertex model 
to a five- vertex model. 

Recently, Verstraete and Cirac have proposed the use 
of continuous matrix product states (cMPS) 14 ) as vari- 
ational states for one-dimensional continuum models in- 
cluding the Lieb-Liniger model. 15 ) It is a natural ex- 
tension of the lattice MPS. The Lieb-Liniger (or the 
quantum nonlinear Schrodinger) model is a model of a 
nonrelativistic quantum field theory and describes the 
one-dimensional Bose gas with the two-body interaction 
V(x,y) — kS(x — y). Theoretical studies of the Lieb- 
Liniger model, such as the exact calculation of correlation 
functions, 16 ) have received considerable attention owing 
to the recent experimental realization of the trapped one- 
dimensional Bose gas. 17,18 ' This model also has a long 
history as one of the exactly solvable models and has 
several different lattice regularizations. While all of these 
models describe the Bose gas at low densities, some of 
them are non-integrable such as the Bose-Hubbard model 
obtained by a natural discretization. 19 ) Integrable lattice 
discretizations can be constructed from the i?-matrix, 
which ensures the integrability through the Yang-Baxter 
relation. 

In this letter, following our previous paper 12 ) , we con- 
struct a cMPS of the Bethe state for the Lieb-Liniger 
model by taking the continuum limit of the MPS of 
the Bethe state in the integrable lattice model. It is 
known that this quantum lattice Hamiltonian obtained 
by the discretization becomes quasi-local, and includes 
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long-range interactions that disappear in the continuum 
limit. 20 ' 21 ) One important consequence of this letter is 
that the F-matrices are indispensable for explicit rep- 
resentation of cMPS, while lattice MPS can be obtained 
without them. Simplification due to the F-matrices leads 
to the five-vertex model for the Lieb-Liniger model, as 
for the Heisenberg chain. To clarify the physical mean- 
ing in the explicit form of the cMPS, we introduce a 
"world line" representation of cMPS as an analog of the 
world line in the continuous time loop algorithm. 22 ' We 
also comment on the cMPS defined by Verstraete and 
Cirac. 14 ' In this comparison, we show that the boundary 
operator 2 ' plays an important role in fixing the number 
of particles in the Bcthe state. 

Let us start from the Hamiltonian of the Lieb-Liniger 
model, which is given by 



product: 



H 



[d x ip^ (x)d x ip(x) + Kip^ (x)^ (x)ip(x)ip(x)] dx, 



where n > and we have fixed h = 2m = 1 . The bosonic 
field operators satisfy the canonical commutation rela- 
tion [^(a;), ijy = S(x — y). For the n-particle state 
|#o) = / d Xl --- dx n ^o(xi,...,x n )tP(x i y ■■■tP(x n y\0), 
one can derive the Schrodinger equation from the above 
Hamiltonian: 



d 2 Xj +2n £ 6( Xj 
i<i<i'<n 



x j ,)}*o = E* , (1) 



where "to = ^0(^1, ■ ■ ■ , x n ). We shall now introduce a 
lattice version of this model, where the spatial position 
x £ R in the continuum model is replaced by the site i e 
Z with the lattice spacing a. Let Vi be a physical Hilbert 
space at the ith site spanned by \m) = -^=f(V'!) m |0) 

with m > 0. Here, ipj and ipi are the bosonic creation 
and annihilation operators on Vi, respectively, and they 
satisfy ipi,ipl = 5u>, and [tpi, ipi>] = tpj,ipl = 0. In 

the continuum limit (a — > 0), ipi — > ij){x)^fa. Note that 
ipi and na are dimensionless. 

Following previous studies, 20,21 ' the L-operator at the 
zth site is denned by 

_ 

/napiipi 
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Kaipjpi 
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where pi = ( 1 + ^fifr-ipi J • The matrix elements of 

Cji and pi are operators on Vi. Let us denote the two- 
dimensional auxiliary space by Vj, which is spanned by 
the two orthonormal states | — >) and | <— ). Cji is repre- 
sented as a 2 x 2 matrix in Vj . This L-operator is the same 
as the classical Lax operator on the lattice and satis- 
fies the Yang-Baxter relation: R(X, X')Cji{X)®Cjn{X') = 
Cj'i(X') (g> £ji(X)R(\, A') with the i?-matrix of the Lieb- 
Liniger model. 20 ' The existence of the i?-matrix ensures 
the complete integrability of this model through the 
quantum inverse scattering method. The monodromy 
matrix is constructed as the following ordered matrix 



N . 

T(\j) = l[£ji(\j)=l 

i— 1 N 



A(A,) 
C(Xj) 



B(A,) 
D(A,) 



where T acts on Vj®J7tf, J7t? = ®i =i Vi is the total Hilbert 
space, and N is the total number of sites. The Hamilto- 
nian defined as i-gjlog tr T(A) is known to be nonlo- 
cal (quasi-local). 21 ' An n-particle eigenstate, the Bcthe 
state, is constructed as 



II B ( A »+i 



;)|o>, 
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where |0) is the vacuum of the total Hilbert space Jff, 
i.e., ipi\0) = for any i. Here, the set of variables 



, n) corresponding to the momenta is the 



solution of the Bethe equation. Since the Bcthe state 
is our starting point in this letter, let us skip the Bcthe 
equation and the algebraic relations among A, B, C, and 
D in the algebraic Bethe ansatz. 

An MPS- type representation of \^> a ) itself is eas- 
ily obtained by changing the order of the product of 
Cji. In eq. (2), first we construct B as B(Aj) = (•<— 

Ilt=i an d t ncn calculate the product of 

B(Aj). Conversely, an MPS representation is obtained by 
calculating the product over j before the product over i. 
Following the notation in ref. 12, the MPS representation 
is given by 



l*»> = Tr 



N 



Qnn( £ *( Ai ' A2 '---' A »)i°» 



(3) 



where _£i(Ai, A 2 , A„) = C ni (X n ) <g> • • • ® £h(Ai), 
jff = V n <g> •• • <g> Vi, and Q n = \ ...,->)(<-,...,<- \ 
= I |. Generally, the MPS is defined by the prod- 

uct of matrices depending on local states. For |\l/ ), the 
local matrix for an m-particle state at the ith site can be 
defined as 

C n , m (Ai,A 2 ,...,A„) = (m|A(A 1 ,A 2 ,...,A„)|0). (4) 

The matrices C n , m , and Q n are 2" x 2" matrices acting on 
The recursion relation between C n +i, m and C n , m ', 
which becomes important for the following discussion, is 
easily obtained as 



si ( bn-\-l,mCn,m Cm— lCn.m— 1 



(5) 



with C„ iTO = (m < or m > n) and C ,o = 1, where 
iXia na 



1 



-m, 



= —i\J (m + l)na (l 



m,Ka\ 



(6) 



(7) 



For the configuration with particles at i\ < i 2 < • • • < i n , 
the component of \^ a ) is given by ^ a (x\, . . . , x n ) — 

Tr^ [QnD^CnD^-^Cn ■ ■ ■ D^-^-'CnD"-' 

where xi = na, D n = C nfi , and C n = C„,i. 

Let us comment on a similar concept, which is the 
Suzuki- Trotter (ST) decomposition 13 ' 23 ' in the context 
of a mapping from a c?-dimensional quantum system into 
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Table I. Correspondence between MPS representation of the 
Bethe state and ST decomposition for an TV-site system. 

MPS for ABA ST decomposition 



ZdWBC = ^a(xi,X2, ■ ■ ■ ,X N ) 
particle number n 
B(Aj) in cq. (2) 
C i (X 1 ,X 2 , . . . , X„) in cq. (3) 



Trotter number n 
T R in Z n = Tr(T R )" 
T v in Z n = TrtTy)^ 



a (d + l)-dimensional classical system. The partition 
function of the two-dimensional classical system mapped 
from the one-dimensional quantum system with N sites 
is expressed as Z = Tr c~P H = Y\m. n =>oo Z n after using 
the Trotter formula, where n is the Trotter number and 
Z n = Tt(Tr)™ = Tr(Tv) N with the real-space transfer 
matrix Tr and the virtual-space (quantum) transfer ma- 
trix Ty. 23 -* Ty paves the way for analytical and numeri- 
cal studies: the thermal Bethe ansatz 24 ) and the finite-T 
DMRG. 25 ' 26 ) 

The similarity becomes clear if we consider a two- 
dimensional statistical model in which the Boltzmann 
weights are given by Cji. For the Heisenberg chain, the 
corresponding model is a six-vertex model with domain 
wall boundary conditions (DWBCs). 27 ) On the other 
hand, in the present lattice model, the partition func- 
tion is defined as the coefficient of the fully filled Bethe 
state as follows: Zdwbc = ^a(xi,X2, ■ ■ ■ ,xn)- Now one 
can find a clear similarity between the MPS of the Bethe 
state and the ST decomposition, as shown in Table I. 
The derivation of cq. (3) from eq. (2) also corresponds 
to that in ST decomposition. The difference is that Z n 
originates from the quantum Boltzmann weight e - ^, 
while Zdwbc originates from the Bethe state \^ a }- 

Hereafter, we shall derive a continuous MPS from ^ a 
defined in the discretized real space with the artificial 
lattice constant a. This reminds us of a remarkable ad- 
vance in the Monte Carlo method, namely, the con- 
tinuous (imaginary)-time loop algorithm, 22 -* which com- 
pletely eliminates the systematic error due to the artifi- 
cial discretization with the Trotter number. 

In the continuum limit a — > with finite n, which 
corresponds to the situation of eq. (1), we need to cal- 
culate an infinite power of D n — C n fl to obtain D z ~ % = 
jj(x x )/a ^ therefore, the diagonalization of D n is in- 
dispensable. To diagonalize D ni we introduce the in- 
vcrtible F-matrices F n and change the basis in ^ as 
Cn, m = F- 1C n,mF n and Q„ = F- x Q n F n . The state |* B ) 
is invariant under this global "gauge" transformation. 
Since D\ is diagonal, we can take the initial F-matrix 
as F\ — cr°, where cr° is the 2x2 identity matrix. Af- 
ter some calculation and using a general procedure (see 
ref. 12), we obtain the following recursive definition for 
the F-matrices: 



+i 



F„ 
F T 







E 



-F, 



(i) 



A.-A 



i A;— An+1 

with the 



and F, 



F n 

where F n 

<7 i f nr=i+i ( j « 

identity matrix I n and 2™ x 2™ matrices a° 



(8) 



x 2" 



®r=i j a ° ® <8>r=n-j+2 defined by the 2 x 2 Pauli 
matrices a a . It is highly nontrivial from eq. (5) that the 
recursively defined F n diagonalizes a non-Hermitian ma- 
trix C n ,o and transforms C n . m >o into upper triangular 
matrices. 

Then, we now write down the explicit expression 



C n ,o — F> n — 



i=i 



bn-l+1,0 






h* 



(9) 



which is diagonal. Other matrices are summarized as 
Q„ = Q n and C„,j = £? =1 C n\ with C n i] = ^C nfl J%\ 
These operators satisfy the following algebraic relations: 

C$V n = zjVnCP, CW)CW = Sjk {Xj , \k)C n k ^C^ , (10) 
where Zi = i 2 ^- and Sa> = z^Sa' with 



Si 



— ~i~ — and Sa' — — 

Moreover, one can show C n . m = 



7-1 



m—1 Cj 



Before we move on to cMPS, let us comment on the 
simplification due to the F-matrices, which enables us 
to make the five- vertex model from the six-vertex model 
for the Heisenberg chain. 12 ) For the present model the 
graphical representation of Zdwbc nas infinite nonzero 
vertices because the local Hilbcrt space Vi has infinite de- 
grees of freedom. Nevertheless, by using the F-matrices, 
we can obtain the five- vertex model under the DWBC. 
In this sense, the present model receives more benefits 
from the F-matrices than the Heisenberg chain. 

Let us consider a diagonal matrix V%~ 11 in the con- 
tinuum limit a — > with a; 2 — x\ = («2 — i\)a. Using 

lim a ^ (l - ia^fj " = e^/ 2 , wc find lim a ^ £>^ 2 ~ 41 = 
exp[i(x 2 - x\)K n ], where JC n = -E™=i^ <7 |- Here ' 



can be interpreted as a free propagator. 



On the other hand, for C n , m (m > 0) one can show the 
relation \fm\C n ,m — (C n ,i) m in the continuum limit with 
a finite n. This guarantees the continuity of the many- 
body function, i.e., ^(x,x) = \im y ^ x ^(x,y). Chang- 
ing the basis by using F-matrices, we can show that 
Ci (Ai , A 2 , . . . , A„ ) 1 0) becomes 



J2Cn, m \m) ~ 2?„exp [ni^(x)a] |0), 



(11) 



where TZ n _= lim a ^ F> n 1 C n ^/^/a= -iv^Eti act- 
ing on Jt° is an analogue of the reflection operator in 

By using the path-ordered operator P, we can write 
down |*o> = lini a ^o |*a) as 



(12) 



The form is seemingly a coherent state with fluctuating 
particle number, but Q„ projects onto the state with a 
fixed number of particles. 

Let us now introduce the graphical representation of 
the cMPS. Following the manner of the Heisenberg chain, 
we write the basis of Vi as | t) = |0) and | \) = |1). The 
m-particle state can be disregarded owing to the conti- 
nuity of the many-body function. Unlike in the lattice 
models, we must draw an infinite number of up arrows 
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Fig. 1. (a) n-particlc cMPS, which is the sum of all configurations 
of "world lines" between momentum space kj = Xj and real space 
Xj. Inset: three connections of world lines, (b) Configuration of 
cMPS with vertices that become zero after using .F-matrices. (c) 
Configurations and corresponding weights after using F-matrices 
for a two-particle state. 



for the vacuum in the continuous space. For simplicity, 
however, we shall draw only down arrows and right ar- 
rows, as shown in Fig. 1(a). 

By using the three kinds of connections shown in the 
inset of Fig. 1(a), the down arrows and right arrows are 
connected one by one, as shown in Fig. 1(b). We call them 
"world lines" in the continuous space and discrete mo- 
mentum space, which is analogous to the world lines in 
the continuous time and discrete lattice. 22 ) To calculate 
$(xi,...,D n ), we sum all configurations of world lines. 
The number of configurations is infinite, because we can 
perform a continuous modification of the configurations, 
as shown in Fig. 1(c). However, after using F- matrices, 
one of the three connections, which corresponds to "anni- 
hilation" , becomes zero. Then, we have only two configu- 
rations for the two-particle cMPS, as shown in Fig. 1(c), 
that is, *(xi,x 2 ) oc e ifelXl e ik2X2 + S(h, k 2 ) e lklX2 e ik ' Xl . 
The simplification due to the F-matrices is valuable for 
developing a new numerical method. 

In the graphical representation, the horizontal line cor- 
responds to the plane- wave- type function e lkx . It orig- 
inates from the free propagator e lfCnX except for the 
overall factor, where K. n is interpreted as a momen- 
tum operator. The meaning of the scattering matrix 
S(ki,k2) = Si2(ki,k2) is clarified by crossed world- lines 
in Fig. 1(c). It originates from the algebraic relation 
eq. (10) and lim a _ ) .o Sjji — Sjj>. In short, the permu- 
tation of plane waves gives the scattering matrix as in 
the coordinate Bethe ansatz. Algebraic relations among 
Q n ,K, n , and Tn\ which define lZ n , are important in 
the continuous matrix product ansatz and can be ob- 
tained generally as a natural extension of those in the 
matrix product ansatz. n ) For an extension of the sta- 
tistical model, a corner transfer matrix (CTM) of the 
cMPS is obtained as a mapping from momentum space 
into real space. Since the original CTM 29 ) is interpreted 
as a Lorentz boost, 30,31 ) it is an interesting problem to 
study the field theory and algebraic structure behind the 
CTM of cMPS. 

Finally, we comment on the cMPS defined by 



Verstraete and Cirac ) variational state: 

|$> = Tr[P(cxp[/M (x)+M 1 (x)V' t (x)dx])] |0), 
where M (x) and M\(x) are variational matrices with 
finite dimension d under the assumption of translational 
invariance for M and Mi, respectively. The exact cMPS 
in the present letter shows finite dimensionality d = 2™ 
and the translational invariance Mq(x) = iJC n and 
Mi (a:) = 1Z n . A significant difference is that while |$) 
shows a particle number fluctuation, the exact cMPS \^f) 
has a fixed number of particles owing to the existence 
of Q n - It is historically interesting that the generalized 
variational MPS 2 ) has the boundary operator Q. In 
short, the exact cMPS has the operator Q n fixing the 
total number of particles as the "boundary" condition 
in the auxiliary space Jf. To tackle future problems 
such as a finite-temperature or higher-dimensional 
generalization, the physical meaning of the matrices 
revealed in the letter would become more important. 
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